Relativistic X-Ray Free Electron Lasers in the Quantum Regime 
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We present a nonlinear theory for relativistic X-ray free electron lasers in the quantum regime, 
using a collective Klein-Gordon (KG) equation (for relativistic electrons), which is coupled with 
the Maxwell-Poisson equations for the electromagnetic and electrostatic fields. In our model, an 
intense electromagnetic wave is used as a wiggler which interacts with a relativistic electron beam 
to produce coherent tunable radiation. The KG-Maxwell-Poisson model is used to derive a general 
nonlinear dispersion relation for parametric instabilities in three-space-dimensions, including an 
arbitrarily large amplitude electromagnetic wiggler field. The nonlinear dispersion relation reveals 
the importance of quantum recoil effects and oblique scattering of the radiation that can be tuned 
by varying the beam energy. 

PACS numbers: 52.35.Mw,52.38.Hb,52.40.Db 
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With the recent development of X-ray free-electron 
lasers (FELs) d there are new possibilities to explore 
matter on atomic and single molecule levels. On these 
length scales, of the order of a few Angstroms, quantum 
effects play an important role in the dynamics of the elec- 
trons. Quantum effects have been measured experimen- 
tally both in the degenerate electron gas in metals and 
in warm dense matter 0, and must also be taken into 
account in intense laser-solid density plasma interaction 
experiments [|[ . The theory of the FEL was originally de- 
veloped in the framework of quantum mechanics Q , but 
where Planck's constant cancelled out in the final FEL 
gain formula producing a classical result. It was sub- 
sequently shown that classical theory can be used and 
quantum effects can be neglected if the photon mo- 
mentum recoil is not larger than the beam momentum 
spread @-[|. To overcome the technical limitation on 
the wiggler wavelength for a static magnetic field wig- 
gler, it was suggested that it can be replaced by an elec- 
tromagnetic (EM) wiggler or by a plasma wave wiggler 
[ici| to generate short wavelength radiation. In such a 
situation, it turns out that quantum recoil effects can be 
important. The Klein-Gordon equation (KGE) for a sin- 
gle electron was used to derive a general set of quantum 
mechanical equations for the FEL [Hj |. while a single- 
electron Schrodinger-like equation for the dynamics of 
the FEL was derived using the field theory [12[ ■ Further- 
more, by using a multi-electron FEL Hamiltonian, it was 
shown that quantum effects can lead to the splitting of 
the radiated spectrum into narrow bands for short elec- 
tron bunches [l3l [l4| . Relativistic and collective quan- 
tum effects have been studied for FELs using Wigner 
. fl(| and quantum fluid 17] models. 

In this Letter, we shall use a collective KGE to derive 
a dynamic model for the quantum free electron laser. 
In our model, we assume that the wave function ip rep- 



resents an ensemble of electrons, so that the resulting 
charge and current densities act as sources [l8| for the 
self-consistent electrodynamic vector and scalar poten- 
tials A and cf>, respectively. 

The KGE in the presence of the electrodynamic fields 
reads 



W 2 tp 



c 2 V 2 ip 



m 2 c A ip = 0, 



(1) 



where the energy and momentum operators are W = 
ihd/dt+ecp, and V = — iKV + eA, respectively. Here fi is 
Planck's constant divided by 2ir, e the magnitude of the 
electron charge, and m e the electron rest mass, and c the 
speed of light in vacuum. The electrodynamic potentials 
are obtained self-consistently from the Maxwell equations 

^+ C 2 Vx(VxA) + v|= W l (2) 



and 



■> 9A 1 , 

at eo 



(3) 



where po is the magnetic permeability, £o = 1/Voc 2 
the electric permittivity in vacuum, and pi = eno a 
neutralizing positive charge density due to ions, where 
uq is the equilibrium electron number density. The 
electric charge and current densities of the electrons 
are p e = -e [<ip*Wip + ip(Wip)*] /2m e c 2 , and j e 

i [tp*Vip + ip(Vip)*] /2m e , respectively. They obey the 
continuity equation dp e /dt + V • j e = 0. 

It is convenient to carry out the calculations in the 
beam frame, and then to Lorentz transform the result 
into the laboratory frame. We assume that the beam is 
propagating along the z-axis, in the opposite direction of 
the laser wiggler beam. For the laser wiggler field, we 
consider for simplicity a circularly polarized EM wave of 
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the form A = (1/2)A exp(— iuj t + ik ■ r)+ complex 
conjugate, with A = (x + iy)A , where uj is the laser 
wave frequency and ko = fcoz the wave vector, and x, 
y and z the unit vectors in the x, y and z directions, 
respectively. Due to the circular polarization, the oscil- 
latory parts in the nonlinear term proportional to A 2 in 
the KGE vanish. Hence, our starting point is the nonlin- 
ear dispersion relation in the beam frame, with primed 
variables, where the plasma is at rest. It reads [19| 



Ik', x A' 



KJ 2 £^(n',K') v 

4 7 3 m 2 c 2 Dj.(n/,K') £ (k' ± rD' A (u' ± ,k' ± ) 



= 0, (4) 



where the electron plasma oscillations in the presence of 
the EM field are represented by 



D' L (Q' ,K') = _ ( n') 2 



1A 

h 2 [c 2 (K') 2 - (n') 2 } 



D' A (n',K'). 



Here fi' and K' are the frequency and wave vector of the 
plasma oscillations, respectively, ^a = (1 + a o) 1/ ' 2 * s the 
relativistic gamma factor due to the large amplitude EM 
field, and do = e\A' \/m e c is the normalized amplitude 
of the EM wave (the wiggler parameter) . The dispersion 
relation for the beam oscillations in the presence of a 
large amplitude EM wave is given by D' L (Cl' ',K') = 0. 
The EM sidebands are governed by 



D' A (u' ± ,k' ± ) = c 2 (k' ± ) 2 - (c^) 2 + « ) 2 /7A, (6) 



where 



± ft' and ki = kn ± K', and 



and fc are related through the nonlinear dispersion re- 
lation (cjg) 2 = c 2 (fcp) 2 + (u' pe ) 2 /jA- We also denoted 
D' A (Q,K) = c 2 (K') 2 - (Vl') 2 + (uS pe ) 2 /lA- We have ne- 
glected the two-plasmon decay [2(| , which would give rise 
to terms proportional to |k' ± • A | 2 in Eq. (j4]). 

To move from the beam frame to the laboratory frame, 
the time and space variables are Lorentz transformed 
as i! — 7o(£ — vqz/c 2 ), x' = x, y' — y and z' = 
7o(z — Vq€), where v = v z is the beam velocity, and 
7o = l/ y/l — Vq/c 2 the gamma factor due to the rela- 
tivistic beam speed. The corresponding frequency and 
wavenumber transformations are thus u) 1 = 7q(w — Vok z ), 
k' x = k x , k' y = k y , and k' z = j (k z ~ v uj/c 2 ), They apply 
to the frequency and wave vector pairs (f2, K), (ujq, k ) 
and (w±, k±). The plasma frequency is transformed as 
u)' pe = u; pe ^7,4/7, where 7 = (1 + pl/m 2 e c 2 + af} 1 / 2 is 
the total gamma factor and po = 7 m e Vo the relativis- 
tic electron momentum. Since the components of Ao 
are perpendicular to the beam velocity direction, they 
are not affected by the Lorentz transformation, hence 
Aq = A . We also use the relation 7 a7o = 7, and ob- 
serve that the expressions of the form (w') 2 — c 2 (k') 2 = 
u) 2 — c 2 k 2 are Lorentz invariant. This yields D' L (il', K') = 



7 2 J Di(ft,K), with 



2 2 

h 2 (c 2 K 2 - n 2 ) , 



(7) 



4 7 2 m 2 c 4 



and D' A (u)'±,k' ± ) = D A (LJ±,k±) = c 2 k\ - uj 2 ± +uj 2 J 1 , 
and, similarly, D A (tt',K') = D A (£l,K). In the labora- 
tory frame, Eq. (U) is of the form 



1 - 



Y— 

4 7 3m 2 c 2 D L (fl,K) (k' ± ) 2 D A {u ± ,k ± 



e 2 |k' + x A n | 2 



= 0. (8) 



Using K = K x 3t + K y y + K z z, we have K 2 = K 2 + K\ 
with K\ = Kl + K 2 , so that |k' ± x A | 2 = {2 7 2 [fc ± 
K z + (v /c 2 )(u; ±n)} 2 + Kl}\A \ 2 , and (k' ± ) 2 = 7 2 [fc ± 



(5) K z + (v /c 2 )(u ±m 2 + KA 



For the resonant backscattering instability, we have 
\D A (u} + ,k+)\ > \D A (L>-,k-)\. Also, for fl « v K z , v « 
— c, and jaK z > ko, Kj_, we have |k' ± x A | 2 /(£4) 2 w 
2|A | 2 . In this limit, Eq. §E§ is written as 



D L {n,K)D A (LO-,k_) = 



u 2 pe D A (n,K) 
2 7 3 



(9) 



By using s» VqK z , the expression for Di(f7,K) can be 
simplified as 



D L (n,K) = n 2 p ~(n-voK z ) 2 , 



where 



47 6 m 2 



(10) 



(11) 



Equation ([TU]) is valid for huj pe ^ A /j 3 ^ 2 m e c 2 -C 1 (which 
is always fulfilled), and w pe 7A/|vol7 3 ^ 2 "C K z <C 
We|wo|7 3/,2 /?i7A- The latter condition, with |«o| 



c, 

gives 7a/A g7 3/2 < X z < j 3/2 /\clA, where A e = c/w pe 
is the electron skin depth and Ac = Ti/m e c the reduced 
Compton length. We note that Q p contains a combina- 
tion of the collective beam plasma oscillation and quan- 
tum recoil effects, which lead to a splitting of the beam 
mode into one slightly upshifted and one downshifted 
mode. 

The resonant Q and K are obtained by simultaneously 
setting D L (Q,K) = and £) A (w_,k_) = 0. Invok- 
ing the approximation Dl(Q,J£) ~ —(fl — vqK z ) 2 = 
and D A (w-,k_) ~ c 2 (fc - K z ) 2 + c 2 K\ - (cj - ^) 2 = 
0, we obtain Q = vqK z and the resonance condition 
(K z -2-f%k ) 2 + j 2 K 2 _ = 4-f 4 k 2 for uj w ck and v « -c. 
The corresponding resonant wave vector components of 
the radiation field, k_ = ko — K, shown in top panels of 
Fig. [31 form an ellipsoid in wave vector space rotationally 
symmetric around the k z - axis. The resulting radiation 
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frequency w_ — uq — vqK z ~ loq + cK z is strongly up- 
shifted in the parallel direction (K± = 0), where we have 
K z = 47 2 fc /7^ and w_ « 47 2 cfc /7A- The result differs 
by a factor two when compared with the case involving 
a static wiggler [21] . 

Comparing the two terms under the square root in 
Eq. (jlll) . we see that the quantum recoil effect starts 
to be important in the parallel direction (Kj_ = 0, 
K z sa 47 2 fc /7^) when k w fc n where 



^O.crit 



^0,crit : 
^pe7j m e 

8fr 7 5 / 2 



1/2 



(12) 

(corresponding to 



In the classical limit kg 
the Raman regime discussed below), we have Q p — 
: while for ko ^> fco, C rit, the quantum ef- 
fects dominate and we have ft p — TiK 2 ^/2j 3 m e = 
&h^k 2 /~f\m e . An expression analogous to (|12l) can be 
derived for the static wiggler case [21] 

For Q «o^, ~ 47 2 fc /7^, i..iu ^ 
have Da(0, K) 4c 2 i4T z fc in Eq. ©, so that 



and w 2 3> w 2 we 



Di(n,K)£l A (a;_,k_) 



2 K z k 



05. 



(13) 



Setting = vqK z + O p + ir, where the real part F# 
of T is the growth rate, and choosing K z and f2 so that 
D A (oj-,k-) = Di(n,K) = for T = 0, we obtain D L = 
-2ift p r + T 2 and Da = 2iT{uj - v K z - Q p ) + T 2 w 
2ir(wo — vqK z ) w 2iTcK z for 3> fco and t>o ~ — c. 
Hence, inserting the expressions for D^ and Da into Eq. 
([IB"]), we have 



r 2 (2^ p + ir) 



CJ Cfc 2 



(14) 



For |F| ^> fl p , we are in the Compton regime where 
the ponderomotive potential of the laser dominates, with 
the growth rate of the instability given by 



T R = 



^3 (w 2 e cfc ) 1/3 2/3 



7 



(15) 



For this case, the quantum recoil effect is negligible [17| . 
On the other hand, for |T| <C \£l P \, we have an instability 
with the growth rate 



(w 2 pe ck \^ 2 



(16) 



Clearly, since Q p is in the denominator, the quantum 
recoil effect leads to a decrease of the growth rate. Com- 
paring Eqs. (|15[) and (|16l) . we find that the limiting am- 
plitude between the two regimes is given by ao = a cr a, 
where 



3o3 \ 1/2 



27 7 ^ 
8w 2 e cfc 



(17) 



Quantum 




FIG. 1: Schematic picture of different regimes for the FEL 
instability, showing the quantum regime fco > fco.crit, the Ra- 
man regime fco < fco.erit, and the Compton regime ao > a cr it- 
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FIG. 2: The maximum growth rate as a function of ao for 
K± = and 7 = 5, using the full dispersion relation |S} (solid 
curve), and the approximations (|14p . (|15[l and (|16[) (dashed, 
dash-dotted and dotted curves, respectively). The vertical 
dotted line indicates ao = a cr i t = 0.038. 



Equation (I15[) is valid for ao 3> a c nt an d Eq. (j!6[) for 
ao -C a cr i t . In the Raman regime ko -C fe Eq. HID 
gives the growth rate w (w pe cfco/27 3 / 2 7A) 1 / 2 a , while 
in the quantum regime fco 3> fco,cHt, Eq. (fl7j| yields « 
(w 2 e 7 A m e c/167 4 ?ifco) 1 / 2 ao. In Fig. 1 , we have illustrated 
different regimes for the FEL instability, including the 
quantum and Raman regimes for small amplitude wiggler 
fields, and the Compton regime for large amplitudes. The 
transition from the quantum to the Compton regime in 
Fig. 1 corresponds to the quantum FEL parameter Q 
p = pm e cy r I TiK z going from smaller to larger values than 
unity, where p = (aow pe /4cfco)/7r is the classical BNP 
parameter [22J, and 7 r the resonant energy in m e c 2 units. 

For illustrative purposes, we choose a beam density 
n = 2.2 x 10 22 m~ 3 , giving to pe = 8.37 x 10 12 s _1 , 
ao = 0.15, and a wiggler wavelength of Ao = 1 y«m giv- 
ing fc = 2tt/A = 6.36 x lO^m" 1 [l6[. For 7 = 5 one 
has ko >cr it — 1.27 x 10' m -1 > fco, so that the plasma 
oscillation effect dominates over the quantum recoil ef- 
fect. Figure displays the growth rate as a function of 
ao, obtained from the exact dispersion relation {5J and 
from the approximations (I14p . (jTH]> and ([15)) . We note 
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k z (m 1 ) x 10 8 k z (m 1 ) x 10 1 ° 

FIG. 3: Resonant parallel and perpendicular radiation 
wavenumbers fc z _ and k±- for 7 = 5 and 7 = 36 (top panels) 
with the corresponding growth rate P(s _1 ) shown in color. 
Bottom panels show the growth rate as a function of k z ~. 
(The scalings of the vertical axes are enhanced in the top 
panels.) The vertical bars show the location of the largest 
resonant wave numbers ^ — 4fco7 2 /7A- 

that the growth rate obtained from (fT4l agrees very well 
with the one obtained from (JSj) . Since a cr u — 0.038 < ao, 
the ponderomotive force dominates over the plasma and 
quantum oscillations, so that Eq. (|15|) can be used to cal- 
culate the growth rate, giving Tr = 2.6 x 10 12 s _1 and an 
interaction length scale c/Tr w 0.1mm. 

On the other hand, due to the quantum recoil effect, 
the gain can rapidly decrease for higher values of 7. Us- 
ing the same parameters as above but 7 = 36 [161 ] . we 
have kQ tCr it — 1.07 x 10 6 m _1 < /co, so that the quantum 
recoil effect dominates the beam oscillations. Here we 
have a crit = 1.6 ;g> a , so that Eq. ([TBI can be used to 
estimate the growth rate, which gives Tr = 1.5 x 10 11 s _1 
and an interaction length c/Tr = 2 mm. For this case, 
the expression (fT5"]) over-estimates the growth rate to 
Tfl/wpe ~ 3.5 x 10 11 s _1 , giving c/Tr ss 1mm. 

The instability of oblique scattering is shown in Fig. 3 
for resonant radiation wavenumbers fc z _ = ko — K z 
and k±- = —K±, obeying the resonance condition 
(K z - 2 7 2 /c ) 2 + 1qK 2 ± = 4 7o 1 fc 2 derived above. The 
growth rate, deduced from Eq. (|14[) . is almost indepen- 
dent of the radiation wavenumbers for 7 = 5, where quan- 
tum recoil effects are unimportant. For 7 = 36, there is 
a significant decrease of the growth rate for larger radia- 
tion wavenumbers due to quantum recoil effects, primar- 
ily in the parallel direction. For too wide electron beams, 



it could lead to a broad-band radiation emission due to 
oblique scattering, while for narrow electron beams this 
is prevented due to a decrease of the possible interaction 
length in the perpendicular direction. 

Summarizing, we have presented a nonlinear model 
for relativistic quantum X-ray FELs, using a collective 
Klein-Gordon model for relativistic electrons, coupled 
with the Maxwell equations for the EM fields, for an 
arbitrarily large amplitude laser wiggler field. We have 
derived a nonlinear dispersion relation for the amplifi- 
cation of the radiation due to scattering instability in 
three-space dimensions. It is found that quantum recoil 
effects can decrease the growth rate of the resonant insta- 
bility, primarily parallel to the beam direction, increasing 
the interaction length over which the radiation amplifi- 
cation occurs. The present study has assumed that the 
coherence of the relativistic electron beam and its trans- 
verse emittance [23| are unaffected by the quantum ef- 
fects over the scale length of out interest. The quantum 
effect could be important if the thermal de Broglie wave- 
length X t h = fry '2tt / 'm e kBT e is comparable to the inter- 
particle distance n^ 1 ^ 3 [24j . For n = 2.2 x 10 22 m~ 3 this 
happens only for T e < 4K when the beam electrons are 
Fermi-degenerate. At room temperature and above, the 
thermal effects clearly dominate over the quantum degen- 
eracy effects on the beam emittance. On the other hand, 
quantum diffusion due to spontaneous photon emission 
could lead to an increase of the energy spread of the elec- 
tron beam. To estimate the relative energy spread, we 
use the formula [H| A7 2 = (14/15)A c r e 7 4 fcgagF(a )Az, 
where F(ao) ~ 1 for ao < 1, r e ~ 2.8 x 10~ 15 m is the clas- 
sical electron radius, and Az is the interaction distance, 
which can be taken to be 10 e-foldings, Az = 10 x c/Tr. 
For the case 7 = 5 above we obtain the energy spread 
A7 = 6 x 10~ 5 , while for 7 = 36 we have the relatively 
large value A7 = 0.014, which might influence the per- 
formance of the FELs. In conclusion, we stress that our 
results will provide a guideline for designing new experi- 
ments for generating tunable X-ray FELs by using tenu- 
ous relativistic electron beams and intense laser wigglers. 
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